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18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting 7^ 6: Laplace's and Poisson's Equations 



We will now study the Laplace and Poisson equations on a domain (i.e. open connected subset) 
c M". Recall that 



(0.0.1) 

The Laplace equation is 



i=l 



(0.0.2) Au(x) = 0, 

while the Poisson equation is the inhomogeneous equation 



(0.0.3) 



Au{x) = f(x). 



Functions u £ C^{il) verifying (0.0.2) are said to be harmonic. (0.0.2) and (0.0.3) are both second 



order, linear, constant coefficient PDEs. As in our study of the heat equation, we will need to supply 
some kind of boundary conditions to get a well-posed problem. But unlike the heat equation, there 
is no "timelike" variable, so there is no "initial condition" to specifyl 

1. Where does it come from? 



1.1. Basic examples. First example: set dfU = in the heat equation, and (0.0.2) results. These 
solutions are known as steady state solutions. 

Second example: We start with Maxwell equations from electrodynamics. The quantities of 
interest are 

• El = (^Ei(t,x,y, z), E2{t,x,y, z), E3(t,x,y, z)^ is the electric field 

• B = {Bi(t, X, z),B2{t^ X, y, z),Bs{t, x, 2;)) is the magnetic induction 

• 3 = (^Ji{t,x,y, z), J2{t,x,y, z), Js(t,x,y, z)^ is the current density 

• p is the charge density 

Maxwell's equations are 



1.1.1) 
1.1.2) 



dtE-yxB 
dtB + VxE 



0, 



V-E: 

V B 



0. 



Recall that Vx is the curl operator, so that e.g. V x B = {pyB^ - dzB2,dzBi - dxB^,dxB2 - dyB^. 
Let's look for steady-state solutions with dtE = 9fB = 0. Then equation (1.1.2) implies that 

1 



2 
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(1.1.3) VxE = 0, 

SO that by the Poincare lemma, there exists a scalar-valued function (f)(x,y,z) such that 



1.1.4) 



E{x,y,z) = -V(i){x,y,z). 



The function is called an electric potential. Plugging (1.1.4) into the second of (1.1.1), and using 
the identity V • V0 = A0, we deduce that 



(1.1.5) 



A(p{x,y,z) = -p{x,y,z). 



This is exactly the Poisson equation (0.0.3) with inhomogeneous term / = -p. Thus, Poisson's 
equation is at the heart of electrostatics. 

1.2. Connections to complex analysis. Let z = x + iy (where x,y e M.) be a complex number, 
and let f{z) - u{z) + iv{z) be a complex-valued function (where u,v e M). We recall that / is said 
to be differentiable at zq if 



(1.2.1) 



lim 

Z-'-ZO 



/(^)-/(^o) 

Z-Zo 



exists. If the limit exists, we denote it by /'(zq). 

A fundamental result of complex analysis is the following: / is differentiable at Zq = Xq + iy^ ^ 
{xQ,yo) if and only if the real and imaginary parts of / verify the Cauchy-Riemann equations at zq : 



(1.2.2) u^{xo,yo) = Vy{xQ,yo), 

(1.2.3) Uy{xo, yo) = -v^{xo, yo). 
Differentiating ( 1.2.2[ ) and using the symmetry of mixed partial derivatives (we are assuming here 



that u(x,y) and v(x,y) are near (xo,yo)), we have 



A dcf 

Am = Urr + U 



'^yx '^xy 



0, 



(1.2.5) V^^ + Vyy = -Uy^ + U^y = 0. 

Thus, the real and imaginary parts of a complex- differentiable function are harmonicl 



2. Well-posed Problems 

Much like in the case of the heat equation, we are interested in well-posed problems for the Laplace 
and Poisson equations. Recall that well-posed problems are problems that i) have a solution; ii)the 
solutions are unique; and iii)the solution varies continuously with the data. 

Let Q c be a domain with a Lipschitz boundary, and let N denote the unit outward normal 
vector to dQ. We consider the PDE 



(2.0.6) Au{x) = f(x), xen, 

supplemented by some boundary conditions. The following boundary conditions are known to lead 
to well-posed problems: 
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(1) Dirichlet data: specify a function g(x) defined on dQ such that u\snix) = g(x). 

(2) Neumann data: specify a function h{x) defined on dfl such that "^7 j^u(x)\dn(x) = h(x). 

(3) Robin-type data: specify a function h(x) defined on dfl such that V Nu{x)\dn(x) + (xu\sn{x) = 
h{x), where co is a constant. 

(4) Mixed conditions: for example, we can divide dQ into two disjoint pieces dfl = Sd u Sn, 
where is relatively open in dfl, and specify a function g{x) defined on Sn and a function 
h{x) defined on such that u\s^(x) = g(x), V^m|sjv(^) = h(x). 

(5) Conditions at infinity: When fl = M", we can specify asymptotic conditions on u{x) as 
1x1 ^ oo. We will return to this kind of condition later in the course. 



3. Uniqueness via the Energy Method 



In this section, we address the question of uniqueness for solutions to the equation (0.0.3), sup- 



plemented by suitable boundary conditions. As in the case of the heat equation, we are able to 
provide a simple proof based on the energy method. 

Theorem 3.1. Let fl c be a smooth, bounded domain. Then under Dirichlet, Robin, or mixed 
boundary conditions, there is at most one solution of regularity u e C'^{fl) n C^{fl) to the Poisson 
equation (0.0.3). 



In the case of Neumann conditions, any two solutions can differ by at most a constant. 



Proof. If u and v are two solutions to (0.0.3) with the same boundary data, then we can subtract 



them (aren't linear PDEs nice?!...) to get a solution w'^= u-v to the Poisson equation with data: 



(3.0.7) 



Aw = 0. 



Let's perform the usual trick of multiplying (3.0.7) by w, integrating over fl, and integrating by 
parts via the divergence theorem: 

(3.0.8) 0= / wAwd'^x= I wV ■Vwd'^x = - f |Vwpc?"x+ / wV<.wda. 

Jn Jn Jn Jan 



In the case of Dirichlet data, w\qq, = 0, so the last term in (3.0.8) vanishes. Thus, in the Dirichlet 
case, we have that 



(3.0.9) 



/ 

Jn 



vwr = 0. 



Thus, \/w - in fl, and so w is constant in fl. Since w is on dfl, we have that w = in fl, which 
shows that u = v in fl. 

Similarly, in the Robin case 



(3.0.10) 

which implies that 



/ w\I<jwda = -(X I w'^da<0, 
Jan Jan 



(3.0.11) f |Vwp = 0, 

Jn 

and we can argue as before conclude that w = in f2. 
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Now in the Neumann case, we have that V jsiwlgn = 0, and we can argue as above to conclude that 
w is constant in Q. But now we can't say anything about the constant, so the best we can conclude 
is that u = V + constant in Q. 

□ 

4. Mean value properties 

Harmonic functions u have some amazing properties. Some of the most important ones are 
captured in the following theorem, which shows that the pointwise values of u can be determined 
by its average over solid balls or their boundaries. 

Theorem 4.1 (Mean value properties). Let u{x) be harmonic in the domain Q c M"-, and let 
Br{x) o Q be a ball of radius R centered at the point x. Then the following mean value formulas 
hold: 



n r 

(4.0.12a) u{x)= u{y)d^y, 

u;„i?" Jbr{x) 

(4.0.12b) u{x)= \ , f u(a)da, 

where ujn is the area of dBi(0) c M", that is, the area of the boundary of the unit ball in M". 

Proof. Let's address the n = 2 case only; the proof is similar for other values of n. Let's also assume 
that X is the origin; as we will see, we will be able to treat the case of general x by reducing it to 
the origin. We will work with polar coordinates {r,6) on M?. For a ball of radius r, we have that 
the measure da corresponding to dBr{0) is da = rdO. Note also that along 9i?r-(0), we have that 
drU = Vm ■ N = \/ j^u, where N{a) is the unit normal to dBr{0). For any < r < R, we define 

(4.0.13) ,(r) ¥ ^l^f^^^^^ uW) = l^f'^ ru(r, 0) - n(r, 0) iO. 

We now note that since u is continuous at 0, we have that 
(4.0.14) u(0) = lim g{r). 

Thus, we would obtain (^TUTT^BJ in the case x = if we could show that g'{r) = 0. Let's now show 
this. To this end, we calculate that 

(4.0.15) g'(r) = ^-l^Mr, 0) dO - ^L^" v^,.(r. 0) M - V^,„«(.) da. 

By the divergence theorem, this last term is equal to 

(4.0.16) Lr Au{y)d^y. 

In Jbi{o) 

But Au = since u is harmonic, so we have shown that 



(4.0.17) 

and we have shown (^TU7T2B} 



g'{r) = 0, 



for X = 0. 



To prove (^ .0.12a) we use polar coordinate integration and ( ^. 0.12b)" (in the case x =0) to obtain 
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(4.0.18) M(0)i?V2= f^ru(0)dr = — ruir.O) dO dr = — f 

Jo 2n Jo Je=o 2n Jb 



u{y)d^y. 



We have now shown (4.0.12a) and (4.0.12b) when x = 0. 



To obtain the corresponding formulas for non-zero x, define v{y) u(x + y), and note that 
Ayv{y) = (Ayu)(x + y) = 0. Therefore, using what we have aheady shown, 



(4.0.19) u(x) = v(0) 



2 r 2 r 

v{y)d^y = — — / u(x + y)d^y = — — / u(y)d^y 

UJ2ri JBr{0) UJoK'^ Jb„(x) 



(^2^^ JBr{x) 



which imphes (4.0.12a) for general x . We can similarly obtain ( |4. 0.12b ) for general x. 



□ 



5. Maximum Principle 



Let's now discuss another amazing property verified by harmonic functions. The property, known 
as the strong maximum principle, says that most harmonic functions achieve their maximums and 
minimums only on the interior of Q. The only exceptions are the constant functions. 

Theorem 5.1 (Strong Maximum Principle). LetQ c M" be a domain, and assume thatu e C{Q) 
verifies the mean value property (4.0.12a). Then if u achieves its max or min at a point p e Q, then 



u is constant on Vt. Therefore, if Vt is hounded and u e C(fi) is not constant, then for every x 6 fi, 
we have 



(5.0.20) 



u(x) < maxu(y), 



u(x) > m.mu{y). 



Proof. We give the argument for the "min" in the case n = 2. Suppose that u achieves its min at a 
point p e Q, and that u{p) - m. Let B{p) c f2 be any ball centered at p, and let z be any point in 
B{p). Choose a small ball Br{z) of radius r centered z with Br{z) c B{p). 
Note that by the definition of a min, we have that 



(5.0.21) 



u(z) > m. 



Using the assumption that the mean value property (4.0.12a) holds, we conclude that 
(5.0.22) 



m - 



1^1 ,1 f uiy) d^y = , ^1 ,1 I [ u{y)d^y+ f u(y)d'^y\ 
\B{p)\Jb(p) ^ \B{p)\\jB.iz) ^ JB\Briz) n 

^ \^B^{z)\u{z) + u(y) d^y] > ^ {mz)\u{z) + m(|5(p)| - \B^{z)\)} 



Rearranging inequality (5.0.22), we conclude that 



(5.0.23) 



u{z) < m. 



Combining (5.0.21) and (5.0.23), we conclude that 



(5.0.24) 



u{x) = m 
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holds for all points x e B(p). Therefore, u is locally constant at any point where it achieves its min. 
Since fl is open and connected, we conclude that u{x) = m for all x e f2. 

□ 

The next corollary will allow us to compare the size of two solutions to Poisson's equation if we 
have information about the size of the source terms and about the values of the solutions on dVl. 



The proof is based on Theorem 5.1 



Corollary 5.0.1. Let c M" he a hounded domain and let f e C(Q). Then the PDF 
(5 25) I An = 0, x € n, 

has at most one solution Uf e C^{Q) nC{Q). Furthermore, if Uf and Ug are the solutions corre- 
sponding to the data f,g^ C(^), then 

(1) ( Comparison Principle) If f > g on dQ and f g, then 

Uf > Ug in Q. 

(2) (Stability Estimate) For any x eQ, we have that 

\uf{x) - Ug(x)\ < max \f(y) - g{y)\. 

Proof. We first prove the Comparison Principle. Let w = uj - Ug. Then by subtracting the PDEs, 
we see that w solves 



(5.0.26) 



Aw = 0, X € 1], 

u(x) = /(x) - (7(x) > 0, X € (9fi, 



Since w is harmonic, since /(x) - g{x) > on and since f g, Theorem 5.1 implies that w is 
not constant and that for every x e Q, we have 

(5.0.27) w{x) > max f{y) - g{y) >0. 

yedfl 

This proves the Comparison Principle. 

For the Stability Estimate, we perform a similar argument for both ±w, which leads to the 
estimates 



(5.0.28) 
(5.0.29) 



w{x) > -max\f{y)-g(y)l 
-w{x) > -max\f{y)- g{y)\. 



Combining (5.0.28) and (5.0.29), we deduce the Stability Estimate. 

The "at most one" statement of the corollary now follows directly from applying the Stability 
Estimate to w in the case f = g- □ 
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